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LIST  OF  SYMBOLS 


" complex  shading  of  the  (m,n)  hydrophone 
tx  - cos  (ax)  - cos  (a^) 

Ty  " cos  (ay)  - cos  (ayR) 


cos  (ax),  cos  (ay) 


cos  cos 


yR' 


direction  cosines  of  radius  observation 
vector 

“ direction  cosines  of  maximum  response 
axis 


propagation  constant 
velocity  of  propagation 


1 
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alx  + aly 


a2x  + a2y  + 2 alx  aly 


“lx’  aly 


expansion  coefficients  for 
i TTNdTyf  | from  Figure  2 


— nMdTxf 1 and 


cos  a. 


salinity  in  parts  per  thousand 


2.34  x 10 


-6 


fT  (KHZ)  - 21.9  x 10  t6  ‘ 152°/(T  + 273)J 


temperature  in  °C 


iv 
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I.  INTRODUCTION 


In  underwater  acoustic  systems,  signal  distortion 
is  often  one  of  the  most  serious  factors  limiting  system 
performance.  Using  the  ambiguity  function  as  a measure  of 
system  performance,  the  following  paper  describes  the  changes 
that  result  in  this  measure  as  signal  distortion  is  included. 
The  sources  of  distortion  considered  are:  (a)  the  sonifying 
and  receiving  arrays,  (b)  the  dispersion  of  the  water,  and 
(c)  the  effect  of  the  inhomogeneity  of  the  water  due  to  the 
small  micro fluctuations  in  temperature. 

Since  both  deterministic  and  random  kinds  of  distor- 
tion are  included,  two  approaches  to  describe  the  effect  on 
system  performance  are  considered.  For  deterministic  distor- 
tion, types  (a)  and  (b)  above,  the  distortion  is  described 
in  terms  of  an  equivalent  transfer  function.  By  relating 
this  transfer  function  to  the  classical, two-dimensional 
ambiguity  function  the  effect  on  system  performance  is  derived. 
For  the  random  distortion  case,  type  (c) , the  measure  of 
system  performance  used  is  the  expected  value  of  the  magnitude 
squared  of  a random  ambiguity  function. 
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II.  DETERMINISTIC  DISTORTION 


Two  sources  of  deterministic  distortion  that  an 
underwater  acoustic  signal  experiences  are  considered  in  this 
section.  The  first  is  that  due  to  the  array;  the  second  is 
that  due  to  dispersion  of  the  medium. 

By  describing  the  frequency  dependence  of  the  array 
amplitude  pattern  and  the  water  dispersion  in  a Taylor  series 
about  the  signal  frequency,  the  effect  of  the  equivalent 
transfer  function  on  system  performance  is  easily  noted. 


Arra^ 


In  both  underwater  acoustic  as  well  as  electro- 


magnetic systems  one  of  the  two  linear  beamforming  beam 
stearing  techniques  is  often  used.  The  first  technique 
utilizes  true  time  delays  for  each  hydrophone,  the  second 
utilizes  phase  shifters.  To  describe  the  frequency  dependence 
of  the  Fraunhofer  region  pattern,  a regular  spaced  planar 
array  is  assumed;  the  extension  to  other  cases  is  straight 
forward. 

\ 

Assuming  the  array  geometry  of  Figure  1 it  is  well 
known  that  the  amplitude  pattern  is  given  by: 
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FIGURE  I ARRAY  DIAGRAM 


From  this  expression  a frequency  dependence  exists 

unless  t - t « o.  Therefore  in  general  one  can ‘express 
x y 

this  dependence  as: 


(f,  V V " E It  cA  (fo’  V V (f"fo)n 

n=o 


For  a dense  uniformly  shaded  linear  array: 


A <£.  O ■ sine  jc  M d Tx  fo  1 + T“  al  ^f_fo^  + TT2  a2  + 


The  coefficients  a^,  a 2 as  a function  of  (i  rr  M d tx  fQ)  are 
plotted  as  Figure  2. 

The  extension  to  the  rectangular  planar  array  is  easily  shown 


to  be: 


r 1 r . 1 

A (f,  Tx,  Ty)  - Sine  M d Tx  f j sine  N d Ty  f j 


which  can  be  expanded  as 


A (f , tx,  Ty)  - sine  £ M 


d 7x  foj  sinc[|N  d Ty  fo]  I1 


+ i-  «,  (£-£  ) + L (£-£„)2  + ...  I ' 


The  usual  phase  shift  implementation  assumes  a phase 
shift  at  each  transducer  output  such  that  at  a frequency,  say 
fQ,  a constant  phase  pressure  wave  in  the  desired  direction 
is  a plane.  For  this  case 


M-l  N-l 


A (f,  Tx,  Ty).  - Y,  t e1  [m(£  cos  ax  - fo  iOS  axR> 


m=o  n=o 


- n (f  cos  a - fQ  cos  a R)] 


Using  the  above  assumptions,  it  easily  follows  that  for  the 
linear  array: 


A (f,  tx)  = sine  — M d (f  cos  ax  - fQ  cos  axR) 


A (f,  Tx)  - sine,  [i  M d Tx  f0]  1 + f.  ax  <f-fQ) 


+ ^ (f-fQ)^  + ... 

o 


The  extension  to  the  rectangular  planary  array  yields: 

A (f,  tx,  t ) - sine  ^ M d tx  fQ  sine  j | N d t f ] / 1 
J J L y J j 

+(fj*l  (f"fo>  + * (*J  *2  (f‘fo)2  | 
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Therefore,  it  is  seen  that  unless  the  maximum  response  axis 
as  well  as  the  observation  vector  are  normal  to  the  array 
(a  p * a ~ 0 for  the  linear  array  and  a D *=  a„D  “ a„  « a..  = 0 


xR  ~ ^x  “xR  yR  ~x  ~y 

for  the  planar  array)  this  array  function  is  frequency  depend- 
ent. Although  the  above  expressions  become  indeterminate  at 
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tx  “ 0,  it  is  possible  to  show  that  a-j_  ■ — j-  a^ 
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Dispersion 


A great  deal  of  work  has  been  directed  to  the  charac- 
terization of  the  dispersive  effects  of  sea  water  on  sound.  It 
is  now  generally  concluded  that  the  significant  factors  include 


the  relaxation  effects  of  the  MgSO^  ions  and  viscosity.  Schulkin 


and  Marsh'1’  have  suggested  that  sound  loss  follows: 


|SA  fT  f 
\f^  + f 
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1 - 6.54  x 10 


Expanding  this  expression  as 

H(f>  - H(f0)  *>„•  - f0>"  • 


' , v 
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the  coefficients  b as  a function  of  f are  shown  as  Figures 
3,  4,  and  5 for  ranges  of  1,000,  10,000,  and  50,000  yards. 

These  assume  a temperature  of  10°C,  salinity  of  35  parts  per 
thousand,  and  zero  atmospheres. 

System  Performance 

Using  the  system  model  as  shown  in  Figure  6,  the 
effect  on  system  performance  of  the  deterministic  signal  dis- 
tortion considered  above  can  be  obtained.  By  assuming  the 
Doppler  approximation*,  it  follows  that  the  equivalent  transfer 
function  is  given  by: 

H (f)  = As  (f-«,  rx,  Ty)  Hf  (f-fi)  (f)  Ar  (f,  tx,  Ty) 


rZAn  <f-Vn 


From  Appendix  I,  it  is  seen  that  the  response  is  given  by 


I 12irf_T  r-*  i r 1 

R (t,  .)  - 2 Re/  r . Y.  AS  T2n]  ~TnT  LX  (T>  4>J 


!n  a<n>  r. 


* For  a discussion  of  this  approximation  see  Remley  . 
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FIGURE  6 SYSTEM  MODEL 
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III.  RANDOM  DISTORTION 

In  addition  to  the  above  sources  of  distortion  is 
the  distortion  caused  by  the  inhomogeneity  of  the  water.  There 
are  small  microfluctuations  in  temperature  which  give  rise  to 
random  fluctuations  in  the  velocity  of  sound  from  point  to 
point  in  the  ocean.  If  one' considers  rays  propagating  through 
these  microfluctuations  the  sound  will  not  be  constant  but  will 
fluctuate  about  the  mean  from  ray  to  ray. 

3 

Liebermann  has  measured  the  fluctuations  of  temperature 

at  depths  of  30-60  meters  using  a fast  acting  thermometer 

mounted  on  a submarine.  He  found  that  the  mean  temperature 

fluctuation  was  0.04°C  and  the  correlation  of  these  fluctuations 

Ixl 

as  a function  of  distance  is  given  by  e * a where  a = 60  cm. 

Chernov^  using  Rytov's  Methpd  was  able  to  derive 
certain  covariance  functions  for  the  phase  and  normalized 
amplitude  fluctuations  for  small  random  deviations  in  the 
refractive  index.  Of  particular  interest  here  are  his  deriva- 
tions of  the  covariance  functions  of  both  amplitude  and  phase 
as  a function  of  distance  both  along  the  ray  path  as  well  as 
normal  to  the  ray  path  (see  Figure  7) . 

For  lateral  (along  the  ray  path)  displacement  of  aL 
he  shows,  using  Liebermann' s results  as  discussed  above,  that 
the  covariance  for  both  amplitude  and  phase  is  given  by: 
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For  the  case  of  transverse  displacements  of  i (normal 
to  the  ray  path)  he  shows  that  the  covariance  for  both  amplitude 


and  phase  is  given  by: 
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The  correlation  between  amplitude  and  phase  is  given 


o - kfi2  loe  ■ 4L  ^ 
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It  is  interesting  to  note  that  if  the  correlation 
intervals  are  defined  as  the  values  of  AL  and  t for  which  the 
covariance  functions  reach  0.5  of  their  maximum  value  we  obtain 


AL  - 0.5  k a4 
t - 0.835  a 
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Therefore,  note  that  decorrel 
slower  than  that  normal  to  th 


'ollows  that  if  one. 


From  the  above  discussion 


assumes  that  the  ray  path  is  changing  over  the  signal  transmission 
as  would  occur  in  the  case  of  a moving  source  receiver,  the 


received  signal  could  be  represented  as 


As  shown  in  Appendix  II,  if  this  signal  is  matched  filte 
detected,  the  response  is  given  by: 


which  pe 


IV.  CONCLUSION 


For  the  deterministic  distortion  case  it  is  seen  that 
the  effect  on  the  system  response  function  can  be  derived  in 
terms  of  coefficients  of  the  Taylor  Series  expression  of  the 
distortion  and  derivatives  of  the  classical  ambiguity  function. 
The  use  of  this  technique  is  attractive  when  the  series  con- 
verges rapidly.  In  addition  curves  were  derived  to  define 
the  coefficients  for  array  derived  distortion  as  well  as 
anomalous  absorption. 

For  the  random  distortion  case  it  was  shown  that  the 
results  of  the  microfluctuations  in  temperature  could  be 
accounted  for  by  convolving  the  resultant  spectral  density 
of  the  error  signal  with  the  classical  ambiguity  function. 


APPENDIX  I 


It  was  found  that  the  equivalent  transfer  function 
for  the  deterministic  distortion  case  could  be  represented  as 


ion 


it  easily  follows  that  the  response  as  a function  of  doppler 
shift  i and  time  shift  t ■ T *t  is  given  by: 


APPENDIX  II 


In  this  section  it  is  desired  to  compute  the  expected 
value  of  the  magnitude  squared  of  the  response  function  R(t,§) 
in  terms  of  the  classical  ambiguity  function  X(t,$).  If  we 
assume  narrow  band  signals  it  easily  follows  that  the  time 
delay  random  variable  is  equivalent  to  a phase  change 
y (t)  * 2nf  y(t)  which  permits  us  to  write; 


[l+XCtpj  |l+X(t2) 


L1+X(ti> 


Now,  by  assuming  that  the  random  processes  are  stationary 
Gaussian  and  making  use  of  the  Wiener-Kn inchin  theorem  yields 
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